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Let us analyze the following evolution equations 
E[x(t), f(t)] — for the variables x(t) (position) and u(t) 
(velocity) of a relativistic particle of mass M > 



M 



du 
~dt 

dx 



u(t), u{0) = u , x(0) = x , 



(1) 



where xq and uq are the initial conditions, 7 > rep- 
resents the damping coefficient and f(t) is a T-periodic 
driving force Notice that defining the momentum 



P(t) 



Mu(t) 



y/l-u*(t)' 

we can transform Eq. ([1]) into the linear equation 

f = -<*•->«• 

where j3 = j/M, whose solution is given by 



P(t) = P{Q)e~ pt - f dzf(z)e 
Jo 



(t-z) 



(2) 



(3) 



(4) 



Equation ([T]) is invariant under time shift (S : t 1— > t + 
T/2) along with the change x i-> — x, provided (Sf)(t) = 
f (t + T/2) = - f(t). The bi-harmonic force 



f(t) = ei cos(qcut + fa) + e 2 cos(put + fa), 



(5) 



preserves this symmetry if, both, p and g are odd integer 
numbers, so in this case the average velocity 



= lim — 

t->+oo t 



u(t) dr, 



(6) 



is zero. In contrast, if p + q is odd and p and q are 
coprimes, a nonzero average current can appear. For the 
sake of simplicity we will take p = 2 and q = 1 in Eq. §5§ 
• Then the solution to (j4} for the chosen force (O will 
be 



P(t) = P exp(-[3t) 
(2 



V^ 2 + - - 



2 + CJ 2 

cos(2w£ + fa - X2), 



cos(u)t + fa — xi) 
(7) 



with Pp = P(0) + (ei/V^+^)cos(^ - xi) + 
( e 2/\//3 2 + 4w 2 ) cos((/> 2 - X2), Xi = arctan(w//3), and 
X2 = arctan (2u>/(3). From ©, one obtains 



tt(t) 



fe=0 



( -l)^(l/2) fc 



2k+l 



(8) 



where (l/2) fe = (l/2)(l/2 + 1) • • • (1/2 + jfe - 1). From 
© and ([5]) it follows that the time-average velocity, v, 
cannot be expressed as a function of the odd moments 
of f(t), unless P(t) is proportional to f(t). Indeed, it is 
only in the overdamped case [in which the inertial term 
in ((T|) is neglected] that the evolution equation is given 
by P(t) = — (1/ 0)f{t) and then v do admit an expansion 
in odd moments of fit). 

Moreover, for small amplitudes ei and £2, the leading 
term of the time-average velocity ([8]) reads 



v = Be\e 2 cos(20i - fa + 9 ), 



(9) 



where B = 3/(8M 3 (/3 2 + uj 2 ) + 4w 2 ) and 6 = 
—2\i + X2- This expression is in agreement with the pre- 
diction of our theory. Furthermore, in the limit j3 — > 
we have — 2xi + X2 7r/2, and in the combined limit 
M — > and /3 — > 00, with 7 = const., — 2xi + X2 —> 0. 
One can check that in the former case Eq. (JXJ) is in- 
variant under time reversal (1Z : t t-> —t) provided 
(TZf)(t) = f(-t) = f(t), and therefore O = ir/2 is the 
prediction of our theory. In the latter case, however, it 
is (Kf)(t) = f(-t) = -f(t) that leaves Eq. (JTJ invariant 
and then our theory predicts 6$ = 0. 
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Equations describing the evolution of particles, solitons, or localized structures, driven by a zero- 
average, periodic, external force, and invariant under time reversal and a half-period time shift, 
exhibit a ratchet current when the driving force breaks these symmetries. The bi-harmonic force 
f(t) — ei cos(qujt + 4>i) + £2 cos(pujt + fa) does it for almost any choice of fa and fa, provided p 
and q are two co-prime integers such that p + q is odd. It has been widely observed, in experiments 
in semiconductors, in Josephson-junctions, photonic crystals, etc., as well as in simulations, that 
the ratchet current induced by this force has the shape v oc e\e\ cos(pfa — qfa + 9o) for small 
amplitudes, where 8o depends on the damping (8q = tt/2 if there is no damping, and 8o = for 
overdamped systems). We rigorously prove that this precise shape can be obtained solely from the 
broken symmetries of the system and is independent of the details of the equation describing the 
system. 
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Ratchet-like transport phenomena, where a 
net motion of particles or solitons is induced 
by zero-average forces, can be observed in many 
physical systems. Such is, for instance, the dc 
current in semiconductors [H-0], the net mo- 
tion of fluxons in long Josephson Junctions (JJs) 
0,0], of vortices in superconductors 0], of cold 
atoms in optical lattices M, \8 1 , or the rectifica- 
tion of Brownian motion [9r|ll| . In some of these 
systems, the ratchet-like motion is induced by 



means of spatial asymmetries [12l . |13|. In the 



others the transport can also appear if some tem- 
poral symmetries are broken by time-dependent 
forces e.g. 13 -Hi]. This latter case has two ad- 
vantages: It is generally easier to analyze theo- 
retically, and it is more amenable to experimen- 
tal observation e.g. in semiconductors |2j, in JJs 
[H, [B[ or in optical lattices 
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A big deal of simulations and experiments have 
show @, i, H, 0, S, [HI, H-S71 that in many differ- 
ent systems the behavior of the ratchet velocity 
v driven by the T-periodic bi-harmonic force 

f(t) = ei cos(qut + 0i) + e 2 cos(put + 2 ), (1) 

where T = 2ti/uj, (pi and </> 2 are the phases, p 
and q are co-primes with p + q odd, and the 
amplitudes €\ and e 2 are small, is given by the 
expression 



where B and 9$ depend on the parameters of 
the model and on u but neither on the ampli 
tudes nor on the phases 0, @, [13, H H H [26 



It has also been shown for specific systems that 
nondissipative dynamics have 9 = n/2 M, 15], 
whereas overdamped ones 

have# = @,EI,E3- 
The aim of this paper is to show that symmetry 
considerations alone are enough to predict the 
behavior (J5J). This is a strong result because it 
is valid for any equation that describes the sys- 
tem, no matter the type of nonlinear terms it 
may contain, as long as it shows invariance un- 
der certain symmetry transformations — which 
will state precisely below. 

Attempts at determining the shape of the cur- 
rent ([2]) can be found even in the pioneering 
works [2 
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Be^el cos(j>0i 



+ #o) 



(2) 



aimed at developing a sensitive 
method of measuring deviations from Ohm's law. 
Their analysis, however, relies on an expansion 
of v in odd moments of f(t), justified by the 
adiabatic response of the system to an applied 
field (see also [28| ) . While it cannot be ruled out 
that such an expansion holds for some systems, 
or in this adiabatic limit, it is certainly not valid 
in general. In fact, if one applies that expan- 
sion to related dissipative systems, like those of 
Refs. @, |, 0, 0, EE 0, 0, M, E}, the value 
#o = is always obtained, whereas 6q ^ in 
general — it can even be #o = vr/2 when dissi- 



2 



pation vanishes. We illustrate this fact by ana- 
lyzing in the Supplementary Material an exactly 
solvable example. There one can readily see that 
the moment expansion is in general an incorrect 
assumption; only in the overdamped or the adia- 
batic limits this expansion becomes correct, but 
we do not know of any proof that this holds for 
systems other than this specific example. 

Let E[x(t), f(t)] = denote a functional equa- 
tion (which can represent an ordinary or partial 
differential equation, an integral equation, etc.) 
describing the evolution of a particle, soliton, 
or localized structure whose position is given by 
x(t), under the driving of a zero-average, exter- 
nal, periodic force fit) = f(t + T), T > 0. One 
such system is said to have ratchet-like behavior 
if the average velocity, defined as [l3[ 

x{t) 



v = lim - 

t^+oo t 



x(t) dr 



lim 

t— >+oo 



t 



(3) 



Si, 



is 



independent of the initial conditions 
nonzero. Consider two temporal transforma- 
tions: time reversal (1Z : t \-> —t) and time 
shift (iS : t H- t + T/2), and suppose that their 
action on the force f(t) is given by 

= /("*) = /(*), (4) 
(Sf)(t) = f{t + T/2) = -f{t). (5) 

Suppose further that any of these transforma- 
tions — with the appropriate transformation of 
x(t) — leaves E[x(t),f(t)] = invariant. Non- 
dissipative systems provide typical examples of 
this kind of behavior such as the equation of mo- 
tions of cold atoms in optical lattices 0] , the dy- 
namic of a particle in a symmetric potential [15] , 
and the soliton ratchets in the extended systems 
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For these kind of systems and forces which sat- 
isfy either (j3J) or (E]) — or both — there can be no 
ratchet effect because any of the two transfor- 
mations changes the sign of v (7Z because time 
goes backwards, and S, because / changes sign). 
As a matter of fact, this is a nice illustration of 
Curie's principle [301 ] . 

In some other cases, time reversal leaves the 
equation invariant provided the force transforms 
itself as 



instead of fll]). The most prominent exam- 
ples of this are equations describing overdamped 
systems such as the vortex motion in JJs jH[, 
the overdamped brownian motion jl6j], and the 
ratchet dynamics of breathers in the discrete 
Schrodinger equation 27 . 

Again, and for the same reason, no ratchet ef- 
fect can appear if the driving force fulfills (E]). In 
this case, however, in general breaking fl^]) and 
(EJ) is not enough to induce a ratchet current, 



some additive noise is necessary as well [21 



Whichever the case, a bi-harmonic force like 
(dD is able to break both @ [or ©] and (J5J) 
and induce a ratchet current. In what follows 
we will prove that, provided a ratchet current is 
produced, the symmetries impose that it be of 
the form j2]). 

Let us begin by noticing that v must be a func- 
tional of f(t), which we can expand as 



;[f] = Vo + j2v n [fi 



(?) 

Vn[f] = (C n (t 1 ,...,t n )f(t 1 )---f(t n )), 



(nf)(t) = f(-t) = -f(t) 



(6) 



where (X) = T~ n J Q dt\ ■ ■ ■ L dt n X and v = 
v[0}. This functional Taylor expansion is a rig- 
orous result of functional analysis valid for a 
very wide class of functionals on Banach spaces 
(see 3l|-33| for details). As v[— f] = —v[f] for 
any force f(t), c 2n (ti, . . . ,t 2n ) = 0, so only odd 
terms appear in the expansion (J7|). On the other 
hand, the functions c n (t 1 , . . . ,t n ) can be taken 
T-periodic in each variable, and can always be 
chosen totally symmetric under any exchange of 
their arguments. Notice in passing that only if 
c n (ti, ...,t n )oz 5(ti-t 2 ) ■ ■ -S^n-i-tn) can v be 
expanded in moments of f(t) — thus the moment 
expansion is only a particular case of ([7]). 

Let us now specialize (J7J) for the bi-harmonic 
force ([1]). First of all, v is not affected by the 
choice of time origin; thus v\T T f\ = v[f], where 
(T T f)(t) =J(t + t) for any r. But f{t + r) = 
ei cosjgt + 0i) + e 2 cos(pt + </> 2 ), with 0i = fa + qr 
and 02 = 02 + pr, so v[f] must depend on the 
phases only through the combination 9 = p<pi — 

<?02 = P4>1 ~ <?02- 

Now we must compute v n [f] for any odd n > 0. 
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By expanding ([T]) in complex exponentials, 



Vn[f] 



ni+n 2 n 3 +n 4 
e 2 



|n|=n, n>0 



x e 



i[(ni-n2)0i+(n 3 -n 4 )02] 



where n = (rii, n 2 , n 3 , n 4 ), |n| = ni+n 2 +n 3 +n 4 , 
and n > denotes a componentwise inequality. 
Besides, because of the symmetry of the func- 
tions c n (ti, ...,t n ), 



A(n) 



where 



<c n (t 1 ,...,t n )e—(' 1 -^)>, 



(9) 



m ™2 n 3 "4 

v = (^7^~?, -g,p7~^7?, -p, • • • , -p). 

The complex number A(n) = i?(n)e^' < ' n - ) , where 
i/}(ni,n 2 ,n 3 ,n4) = -i/)(n 2 , m, n 4 , n 3 ) 



B(ni,n 2 ,n 3 ,ri4) = B(n 2 , n u n 4 , n 3 ). 



(10) 



Let us focus now on the complex exponential 
in (jSJ). We know that i> must be a function of 6*, 
so the only nonzero terms in this sum are those 
satisfying (rii — n2)0i+(n 3 — n 4 )0 2 = k6 for some 
k G Z, i.e. ni — n 2 = fcp and n 4 — n 3 = fcg. But 
that means \k\(p + q) = \n\ — n 2 — n 3 + n 4 | < 
|ni + n 2 + n 3 + n 4 | = n. 

Suppose n < p + q; then A; = 0, which implies 
ni = n 2 and n 3 = n 4 , and therefore n must 
be even. Thus v n [f] must be zero for any odd 
n < p + q, but since there are no even terms 
in (J7|), this means that no term with n < p + q 
contributes to v. 

Suppose now n > p + q; then \k\ > (k = is 
excluded because it would lead to an even n) and 
so there will be nonzero terms in (jSJ) correspond- 

2no + \k\p 2n% + \k\q 2n\ + \k\p 2riA + \k\q 

mg to powers e x e 2 or e x e 2 

The lowest order is n = p + q, and is obtained ei- 
ther when n 2 = n 3 = 0, n\ = p, and n 4 = g (i.e. 
fc = 1), or when ri\ — n 4 = 0, 77-2 = p, and n 3 = q 
(i.e. fc = —1). Because of ffTUl) . the contribution 
of these two terms to (EI) is 



Jp+q 



Be\e\ cos(6 + O ) 



11] 



where B = 2B(p, 0, 0, q) and fl = ^(p, 0, 0, g). 



Let us now assume that the equation 
is invariant when the force satisfies (jlj). 
Then, since (cn(ti, . . . ,t n )f(-ti) ■ ■ ■ f(-t n )) = 
(c n (-t 1 ,...,-t n )f(t l )---f(t n )), and v[Kf] = 
—v[f], hence 



Cn(-h 



~tr 



-c n (h 



(12) 



Applied to © this means that 5(n)e _i ^ (n) = 
-S(n)e^W, i.e. ^(n) = tt/2 for all n. 

On the other hand, if the equation is invari- 
ant when the force satisfies (JS]), then v[—TZf] = 
—v[f], so (recall that n is odd) 



C?i( t±j . . . j t n ) C n (ti, • • • , tn)- 



(13) 



Applied to © this means that £(n)e^ n ) = 
5(n)e^ n ), i.e. ^(n) = for all n. 

What we have just shown is that the mean ve- 
locity, v, of the ratchet current induced by the 
bi-harmonic force ([T]) in an equation which is 
invariant under (|5]) always has the form (j2J) if 
the amplitudes e± and e 2 are small. The coef- 
ficients B and 6q depend on the frequency and 
on the remaining parameters of the system, but 
not in a universal way that can be predicted un- 
der symmetry arguments like these ones. This 
shape for the current has been observed, mostly 
for p = 2 and q = 1, in experimental, numerical, 
and theoretical results in several seemingly unre- 
lated systems 0, 0, EH, 0, MM, Q • For p = 4 
and q = 1, the collective coordinate on soliton 
ratchets developed in 26[ also confirms (T5]). 

If the equation is also invariant under (j3J) 
0, EE)], then 6q = n/2 and we recover the form 
v ~ ef e 2 sin 9, whereas if the equation is invari- 
ant under (jBJ), then 9 = and v ~ e^elcosO, 
in agreement with the vortex motion observed in 
JJs [5j, with the overdamped stochastic dynamic 



of particles studied in [111 . |21| . and with the 



ratchet mobility of breathers in the discrete non- 
linear Schrodinger equation computed for p = 2 
and q = 1 in [27]] . 

Notice that formula ([2]) does not imply that 
B must be nonzero (Curie's principle). It only 
proves, under symmetry arguments, that the 
leading term of v can be of that precise form. It 
might well happen, for some specific equation, 
that B = 0. In this case this analysis shows that 
the leading term must have a dependence on the 



4 



amplitudes through powers higher than p and q. 
It is likely that if this occurs it will be the fin- 
gerprint of a hidden symmetry which, properly 
broken, will restore the result (j2J). 

This analysis provides a direct way to quanti- 
tatively relate the causes and the consequences 
of phenomena through Curie's principle. For in- 
stance, our study can be extended to the so- 
called gating effect, i.e. when the amplitude of 
spatial or field potentials for particles or soli- 
tons, respectively, is modified by a multiplica- 
tive force g(t) = e 2 cos(pojt + 2 ) as well as an 
additive force, f(t) = e\ cos(qut + (bjY with p 
and q coprimes, acts on the system [34U36|. In 
such systems, if both f(t) and g{t) satisfy (jl]) (in 
the non-damped limit) or ([H]) (in the overdamped 
limit), or f(t) fulfills (j3J), a ratchet transport can- 
not be induced. A similar procedure shows that, 
when these symmetries are broken, the average 
ratchet velocity is also given by Eq. (j2J), where 
#0 = or # = 7r/2 in the non-damped or over- 
damped limits, respectively (37|. 
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